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Abstract 

The electromagnetic scattering of a spin-0 charged particle off a fixed 
center is calculated in first-order quantum perturbation theory. 

This implies evaluating the square of a 'Dirac delta-function,' an op- 
eration that is not defined in Schwartz distribution theory, and which in 
elementary text-books is dealt with according to 'Fermi's golden rule.' 

In this paper these conventional calculations are carefully reviewed, and 
their crucial parts reformulated in a Colombeau algebra — in which the 
product of distributions is mathematically well defined. 

The conclusions are: (1) The Dirac delta-function insuring energy con- 
servation in first order perturbation theory belongs to a particular subset of 
representatives of the Schwartz distribution defined by the Dirac measure. 
These particular representatives have a well-defined square, and lead to a 
physically meaningful result in agreement with the data. (2) A truly consis- 
tent mathematical interpretation of these representatives is provided by their 
redefinition as Colombeau generalized functions. This implies that their 
square, and therefore the quantum mechanical rule leading from amplitudes 
to probabilities, is rigorously defined. 



1 Introduction 



Quantum theory and quantum field theory are plagued by fundamental conceptual 
problems which originate in part from the mathematical difficulties of properly 
interpreting and dealing with the '5-functions' which were first introduced in these 
theories by Dirac in 1927, i.e., [1, p. 624], as a development of Lanczos's field 
theory [i2|. 

In this paper we calculate in first-order time-dependent perturbation theory 
the cross-section of what is possibly the most simple interaction between two 
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charged particles: The Rutherford scattering of a spin-0 particle of mass m from 
a fixed Coulomb potential, i.e., the electric field of a nucleus of charge Ze and 
mass M oo. The scattering center being fixed implies that the modulus of the 
incoming momentum pi of the particle is equal to that of its outgoing momentum 
Pf, i.e., \pi\ = \pf\. Its final energy Ef is thus equal to its initial energy Ei, and 
the only effect of the interaction is to change the direction of the incident particle. 

All calculations are made in full detail, and at the most elementary level of 
quantum field theory. 

The objective is to clearly show how (5-functions and their squares arise in that 
problem, and how these 6- and 5^-functions are evaluated in conventional theory, 
as well as in Colombeau's theory of nonlinear generalized functions. 

For an introduction to Colombeau generalized functions see [3], for a compre- 
hensive textbook dUl, and for a primer |[5l or [|6ll. 

2 Definition and normalization of states 

The amplitude (i.e., the state, or wave-function) of a free spin-0 particle is defined 
as 

(j) := uexp{ip ■ X — iEt), (2.1) 

where u E C The probability current density associated to this amplitude is the 
four-vector 

J^ = t[(l)*id^(l))-id^(l)*)(iyy (2.2) 
so that the probability density, given by the time component of that four- vector, is 

^ = ^(^^*{^t<P)-{^t<Pn<p)■ (2.3) 

The state (p is assumed to be an element of a Hilbert space in which the Lorentz 
invariant scalar product is defined as 

((0ill02)) ^= ^ - (2.4) 

The condition normalizing the amplitude (12.11) can therefore be written 

V = j£dV = {{m) = l, (2.5) 
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where V C M'^ is a large 'box' of volume V in which the scattering experiment is 
supposed to take place. Substituting from (|2.1I) this yields 

r = i jjj £x u*u{-iE - iE) = 2u*uEV = 1, (2.6) 

which implies that the normalized amplitude is 

1 



exp(ip ■ X — lEt + i^po), (2.7) 



^J2EV 

where the constant phase v^o will be set to zero in the following. 



3 Calculation of the transition amplitude 

In first order perturbation theory the transition matrix element from a state </)j to a 
state 0/ is given by 

Sfi = 2e Iff d^x [ dt c^}A^i, (3.1) 

where A is the operator corresponding to the perturbation. In our case A is simply 
the Coulomb potential of the nucleus, i.e., in Heaviside units, 

A = (3.2) 



Ze^ 1 ///" di?'x ( \ f^°° ( \ 



Thus, substituting the normalized initial and final state amplitudes, 

ze^ 1 rrr d^x - ^ /-^^^ 

(3.3) 

The spatial integral is the Fourier transform of the Coulomb potentiallll well known 
to be 

d^X /./^ ^ - ^\ 4:71 

-j^exp{^i{Pi- Pf) ■ xj = (3.4) 

where q = pj — Pi is the transferred momentum. The integral over the time coordi- 
nate is not well defined because its integrand is non zero at infinity. Nevertheless, 



'Our conventions are such that the Fourier transform is f{p) = exp(ipx) f{x) dx and its 
inverse f{x) = (l/27r) J^^ exp{-ixp) f{p) dp. 
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in conventional calculations, it is interpreted as yielding the (5-functior|^ 

f dt exp(i{Ef - Ei)t \ = 2tx 5{Ef - Ei), (3.5) 

J — oo 

which expresses energy conservation between the initial and final states in a static 
potential. Thus, with q = \q\, the transition amplitude is 



Remark 1 (Conventional interpretation of the ^-function) Eq. (|3.5I) is mathe- 
matically ill-defined because exp(ix) is undefined at x = ±00. The conventional 
interpretation in quantum theory assumes that there is an implicit cut-off such that 

I r+T/2 ^ .+T/2 

6{x) := — lim / dp exp(±ixp) = — lim / dp cos(xp), (3.7) 
is real-valued and is even in the variable x. 



4 Calculation of the transition probability 

The transition probability dWfi is obtained by multiplying the modulus squared of 
the transition amplitude by the phase-space element, that is by the number of final 
states in momentum interval d'^pj, which is equal to Vd^pf / {2nY. Therefore 

After using (13.61) and making the substitution — > Aixa, where a = 1/137 is the 
electromagnetic coupling constant, this transition probability becomes 

dW,. = 8n^l,^\6m-E^). (4.2) 

in which the modulus squared of the 5-function arises. 



^What is meant here by '^-function' or 'Dirac function' is the loosely-defined object known 
under these names by the physicists, not the mathematically well defined 'Dirac measure' or 
'(5-distribution' of Schwartz distribution theory. 
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The time-integrated total cross section S is defined by the phase-space in- 
tegrated transition probability divided by the incoming flux, that is divided by 
Vi = \vi\/V where Vi = pi/Ei is the incident velocity. Thus 

E := [[[ dWf,-, (4.3) 

and therefore 

S = 8vrZV/// ^^/^mEf-E^). (4.4) 
JJJrs EiVi Ef 

Then, with ViEi = pi, and d^pf = d^lp'jdpf. 

Finally, since by definition E"^ — p^ = m? we ha.\e pfdpf = EfdEf, and thus 

§ = r ^^^JSm - E.)dEf, (4.6) 



where, since 5{Ef —Ei) ^ only for Ej ^ Ei, the integration range was extended 
from [0, -l-oo] to [— oo, +cxd] for convenience. 



5 Calculation of the differential cross section 

Equation (14.61) contains which like 5^ is mathematically undefined if 5 is 
interpreted as a Schwartz distribution. In standard quantum theory this fact is 
generally ignored and 6 is assumed to be a real valued function defined by a 
regularization consistent with Remark[il The conventional method for calculating 
(|4.6I) can then be introduced in several ways. A purely heuristic approach consists 
of blindly generalizing the usual sifting property of the 5-distribution, i.e., 

/+00 
f{Ef - Ei) 6{Ef - E,) dEj = /(O), (5.1) 
-oo 

to the (5^ -function by writing 

/ + 00 
f{Ef-E,)5\Ej-E{)dEf 
■oo 

= (fiEf - Ei)5{Ef - E,))siEf - E,)dEj = f{0)5i0). (5.2) 

J —CO 
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With 5{E) defined according to (13.51) . the diverging quantity 5(0) is then estimated 
by introducing a symmetric cut-off to regularize the integral as in (|3.7I) . i.e., 



'+T/2 



dt exT)(iEt) = — lim 



dt exp{iEt), 



(5.3) 



■T/2 



SO that 



S(0) = lim — . 
^ ^ r^oo27r 



(5.4) 



Applying (|5.2I) and (15.41) to (14.61) . and dividing both sides by T, the time-integrated 
cross-section becomes the differential scattering cross section 



which is known as the Rutherford scattering formula for a spin-0 projectile. This 
result is very well confirmed by experiment as the first term of a perturbation 
series El For this reason, any regularization of the 5-f unction which like (15.31) 
yields (15.41) for 5(0) is a suitable representation of the (5-function for this quantum 
mechanical problem. 

It should also be remarked that the physically measurable quantity, the cross 
section (15.51) . is a rate 'averaged' over a very long time T, not an instantaneous 
rate in a infinitesimal interval dT. Actually, the time T is the 'duration' of the 
measurement. Its product with the the volume V of the 'box' used to normalize 
the amplitude defines a space-time volume VT such that the integration over 
infinite space-time in (13.11) is actually restricted to that finite hyper-volume. 

This heuristic derivation can be improved in various ways, see for example ^ 
p. 101]. Moreover, the same method can be applied to other first order perturbation 
problems for which a general formula, called 'Fermi's golden rule,' can be derived, 
see for example BSl p. 332]. Finally, the method can be further generalized to four- 
dimensions and thus to invariant perturbation theory, e.g., \Jj p. 112], \^ p. 164]. 

6 Justification of the conventional method 

In this section we show that the 5-function regularized according to the standard 
prescription, given in Remark[il is consistent with Eqs. (15.11) and (15.21) . 

''if the projectile were a spin-1/2 particle Eq. i5.5l would be multiplied by a correction factor 
that would reduce to 1 in the non-relativistic limit. 



dn ■ 



da 




(5.5) 
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Let us rewrite (15.21) in the form 

»+oo 



f{x)6\x)dx = 6{0) / f{x)S{x)dx 



+00 



(6.1) 

where, replacing T/2 by 1/e with < e < 1 and making an elementary integration, 

1 sin(x/e) 



the 5-function defined by (15.31) corresponds to 



1 r'^' 

6{x) = lim — / exp{—ipx) dp = lim 

e^O 27r e^O 

Because, [Hoi p. 405], 

+°° 1 sin(a;/e) 



TC X 



dx 



the sequence 



Pe 



-p{-h 

e V e ^ 



TT X 



where 



p{x) 



1 sinfx) 



(6.2) 



(6.3) 



(6.4) 



TT X 

is normalized to 1. However, since has non-compact support, we have to 
carefully verify under which conditions is a (5-sequence. Thus we take f E V, 
i.e., a smooth function with compact support, and calculate 



f{x) 6{x) dx = lim — 
= lim — 



00 
+00 



X 



f{ez) dz. 



(6.5) 



Then, if we try to develop f{ez) in a Taylor series, we are immediately confronted 
with the difficulty that the moments 

"^"^ 1 sm(z) 



TT Z 



2;^" dz = 00, 



Vri > 1. 



(6.6) 



We therefore try another standard method [ | 11. p. 29], which consists of splitting 
the integral into three pieces, i.e.. 



f{x) 6{x) dx = lim — 

e->0 TT 



+ 



'^)f{xf^^dx. (6.7) 



X 



The first and the third pieces tend to zero because f E V implies that there is a 
bound B such that 



lim — 

e^O TT 



±00 



< 



5 lim — 



±00 



sm 



dz 



0. (6.8) 
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The second piece can be computed using the mean value theorem because sin (2;) / z 
is even. Thus, with ^ E [— -^/e, +-\/e], 

n L r-^;w!i£(£Zl),..u r^!!f^,..;,o,. ,6.9, 

Therefore, provided f E 'D,5 defined by (|6.2I) has the sifting property of the Dirac 
distribution, i.e., Eq. (15.11) . thanks to the normalization (|6.3I) . Similarly, since IfTOl 
p. 446] 

pI{x) dx = / — = — ' (6-10) 

J-00 ^ Tre 

the same proof works for the square of b^. Indeed, 

/ /(x)(52(x)rfx = lim— / /(o;)^ — fix = lim — /(O), (6.11) 

implies that the sequence Txtp\{x) also defines a 5-function. Moreover, 5(0) = 
lime^o shown by (|6.2I) . The right-hand side of (|6.1 II) is thus equal to 

that of (|6.1I) . which is therefore proved. 



Consequently, the conventional calculations leading to the Rutherford formula 



(|5.5I) are correct. However, these calculations are entirely based on the special 
properties of the particular representative (16.21) of the (5-function. More precisely, 
for f E V, the proof of (16.11) requires that the (5-sequence p^ must be such that p 
has the two properties 

/•+00 

p{z)dz = l, (6.12) 



and 

" + CXD 



/ + 00 
p\z)dz = p{0), (6.13) 
■00 



and the Rutherford formula is obtained if moreover 



p(0) = -. (6.14) 

TT 



In conclusion, the 5-function defined by (13.51) and regularized according to 



(13.71) is a special representative of the Schwartz ^-distribution having the additional 
properties (16.131) and (|6.14l) . The mathematical object denoted '5' in Eq. (14.61) . 
and named '5-function' by the physicists, is therefore not a Schwartz distribution 
in the sense that its representatives p^ are not any element of the equivalence class 
of representatives corresponding to the functional 5{f) defined by Eq. (15.11) . 
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7 Remarks on the conventional method 



1. The 5-function (16 .21) can be rewritten as 



1 



+00 



5{x) = limpe(x) = lim — / exp{—ipx)p{ep) dp (7.1) 

where instead of being implemented by finite bounds in the integration range, 
the 'sharp' regularizing cut-off is implemented by the Fourier transform of 
p{x), i.e., by the piece- wise continuous function p(z) = 1 for |z| < 1, and 

p{z) = for \z\ > 1. 

Similarly, the conditions (|6.12l) and (|6.14l) on p, i.e., 

/ + 00 
p{z) dz = 1, and p(0) = -, (7.2) 
■oo TT 

can be implemented on its Fourier transform p, i.e., 

/ + 00 
p{z) dz = 2. (7.3) 
-oo 

Thus, any cut-off function p G M satisfying (17.31) as well as (16.131) would 



yield a regularization of (|3.5l) leading to the Rutherford formula, and more 
generally to any formula derived by means of Fermi's golden rule. 

For example, the exponential damping factor p^ = exp(— |z|) satisfies the 
conditions (|7.3I) . It provides a 'smooth' cut-off such that 



1 /"'^°° 1 e 
Pei^) = 7^ dp exp{-ipx) exp{-e\p\) = 5—^, (V.4) 



which is a well-know representative of the (5-function having the property 
(5(0)f = l/(7re). Unfortunately, pf does not satisfy (|6.13l) . As a second 
example, if the regularization is made by means of the Gaussian damping 
factor fP = exp(— 2:^/4), one gets another well-know representative of the 
^-function, i.e., 

Pe = ^ J exp{-ipx) exp(^ — j = -^expi^- — 

(7.5) 

which however does neither satisfy (|6.3I) nor (16.141) because S{0)f = 1/ y/ne. 
It is therefore not easy to find a simple example of a smooth damping function 
satisfying all three conditions (l6.12l[6T3l and l6.14l) . 
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2. In the notation of Sec.[5l where 1/e corresponds to T/2, the conditions (17.31) 
read 

/•+°° /2t\ 

p(0) = l, and / pi-\dt = T. (7.6) 

J —oo ^ / 

Therefore, the physical interpretation of the second condition on p, which 
in the form (17.21) corresponds to the point value p(0) = l/vr, is simply that 
its time integral is equal to T, the effective time over which the scattering 
process is averaged (or measured). 



3. To prove (|6.1I) we assumed that f E V, i.e., that / was smooth and had 
compact support. However, only the facts that / G and f < B where 
actually used in the proof. Let us verify that the function f{E) implied by 
(|4.6I) satisfies these conditions. We have 

f{Ef,Ei) = ^\, (7.7) 
Pi 



where p = ^/E"^ — w? and where, with Q the angle between and pi, 

q\Ej, E,) = \pj -p,\''=pj+ p1 - 2pfp, cos 9. (7.8) 
Thus, for fixed Ei ^ m, i.e., pi ^ 0, the function 

* Pi{p} + pf - "ipfPi cos ey ' 

is C°° and bounded, unless Ej = Ej and = 0. Consequently, excluding 
the latter case (which arises in the limit where there is no interaction) / does 



indeed satisfy the conditions under which (16.11) is valid. 



8 Rigorous method 

There are several way of treating the physical problem addressed in this paper in 
a Colombeau algebra. The most comprehensive approach would be to define all 
objects as elements of a Colombeau algebra right from the beginning, as is done 
for example in [|12 | | . This would provide a formulation having all the advantages 
of the 'wave-packet' descriptions of scattering processes [8^ p. 384]. However, 
for the present problem, it is enough to suppose that the damping function p 
employed to regularize the 5-function (17.11) is a suitable damper insuring that the 
(5- sequence p^ is an element of a Colombeau algebra^ This is sufficient for the 

'^For the definitions of suitable moUifiers and dampers in a Colombeau algebra, and for a 
discussion of their basic properties, see. e.g., llT2l Sec. 1.4]. 
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modulus squared of the (5-function arising when squaring Sji given by (13.61) — 
e.g., when calculating the probability (14. 1|) — to be mathematically well defined 
and meaningful. 

The conditions for p to be a suitable damper are simply that 

peV(R), and p(0) = 1, (8.1) 

which by Fourier transformation are equivalent to the well-known Colombeau 
constraints on the moments of p, i.e.Jfl 

pe5(C), Jdzp{z) = l, and J dz z^'piz) = 0, Vn = l,...,g, (8.2) 

where g G N is as large as we please. 

To these conditions we need to add the special constraints (16.1311041) . i.e., 

j \pWz) dz = p(0), and p(0) = p*(0) = ^, (8.3) 

which insure that the properties of the damper p acting as a 'smooth cut-off are 
equivalent to thoses of the 'sharp cut-off introduced in (15.31) with regards to the 
present problem, i.e., to obtain the Rutherford formula. Note that, since p G M, the 
second condition on the right of (|8.3I) is consistent because by Fourier transform 
p(0) = (27r)^^ J p{p) dp is necessarily real, even though p{x) may be complex for 

To find out the full consequencies of working in the Colombeau setting we 
reconsider (16.11) in which 5{x) G C defined by the representative sequence p^ is 
now subject to the Colombeau constraints (|8.2I) and to the physical constraints 
(|8.3I) . That is, we study 



Um / f{x)\pe\ (x) dx = \im p^{0) / f{x)p^{x)dx, (8.4) 

e^O J e^O J 

where for the sake of generality we suppose that / G OmO^, C), i.e., f{x) G 
and each of its derivatives do not grow faster than a power of x at infinity in M. 
Thus, making the change of variable x = ez, we want to find out under which 
conditions we have 

lim / f{ez)\p\'^{z)dz = Mm p{0) / f{ez)p{z)dz. (8.5) 

e-*0 J e^O J 



^Note that if the damper p e R the molUfier p may nevertheless be complex. 
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We begin with the right-hand side of (18.51) and use Taylor's theorem with 
remainder. We obtain 



limp(O) / piz)f{ez) dz 

: p(0) lim £^ p{z) (/(O) + ezf'iO) + (0) + ...) dz 

:p(0)/(0) + limp(0) / p{z)p^/''+'\^z)dz 
J-oo W + 

:p(0)/(0) + limO(e'+i), (8.6) 



where we used the Colombeau constraints (|8.2I) implying that all moments of 
order n with l<n<g + lare identically zero, and where -(9 e]0, 1[ in the 
remainder. Then we observed that since p E S and /('?+^) g Cm their product 
with z'^^^ is still in S, and therefore the remainder is bounded so that (18.61) is equal 

to /(O) + 0(e''+i). 

Turning to the left-hand side of (|8.5I) and making again a Taylor development 
we get, 

r+oo 

lim / \p\'^{z)f(ez) dz 



lim 



■oo 

+00 



oo 

+ 00 



p\\z) (/(O) + ezfiO) + ^r(O) + ..)dz 

/+00 r+oo 
\pWz) dz + lim(ef'{0) / z\pWz) dz + ... + 0{e''+^) 
■oo ^-"OV J^^ 

p(0)/(0) + lim(er(0) ^ " ^|pp(z) rfz + ... + 0(6"+^)) , (8.7) 



where we used (18.31) to get the last line, and where we did not write the remainder 
explicitly because its form is similar to that in (|8.6I) with p replaced by |pp, so 
that it is again 0(e'^+^). Equation (18.71) cannot be further reduced unless specific 
constraints are put on the moment of |pp. Due to p G 5, the only thing that can 
be said is that these moments are finite. 



Comparing (18.61) and (18.71) we see that equation (|8.4I) can in principle be 
obtained at two different orders of approximation: 



1. Eq. (18.41) is an equality 0(1) a* e — > 0. 

This case corresponds to the conventional perturbation theory calculations 
in which the results are limits such as (15.51) where it does not matter how 
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fast the limit is approached. This is a suitable method if the result is not 
to be used in a subsequent calculation in which it could be multiplied by a 
quantity which might be diverging so that the 0(1) remainder can be safely 
ignored relative to the leading term which is 0(l/e). If that is the case, 
no special constraint on the moments of |pp with n ^ is required, and 
the only constraints on the damper p are those implied by (|8.1I) and (18.31) . 
Expressing them as conditions on its inverse Fourier transform p, they are 



2. Eq. (18.41) is an equality 0(e'^)/or all q eN provided e g]0, 1[. 

This corresponds to the case where (|8.4I) would be an equality in the 
Colombeau algebra, which implies that the difference between the two sides 
of (|8.4I) must be O(e^). Referring to (18.71) this requires that all the moments 
of IpP with n = l,2,...,gare zero. Therefore, the full set of conditions on 
p is 



However, as |pp G M, it is immediately seen that the last condition on 
the second line of (18.91) cannot be satisfied when n is even, unless p is 
identically zero. Consequently, the option "Eq. (18. 4|) is an equality 0(e'^)" 
is not feasible. 



9 Discussion 

There are three major mathematical problems in calculating the cross-section of 
the scattering process considered in this paper: 

1. The interpretation of the time-integral in Eq. (13.31) as a Dirac (5-function 
insuring energy-conservation in the transition amplitude (13.61) : 

2. The interpretation of the square of that (5-function in the transition probability 
(|4.2I) . which derives from the quantum-mechanical rule (14.11) associating a 
probability to the modulus squared of the amplitude of the process; 




(8.8) 




(8.9) 
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3. The integration of that 5^ -function in the differential cross-section (|4.6I) . 



The conventional solution to the first problem is to introduce a sharp cut-off to 
regularize the time-integral in Eq. (13.31) . which therefore (as explained in Remark!!]) 
leads to the particular representation (13.71) of the Dirac measure. This regularization 
is physically well motivated by the notion that a free particle corresponds to a 
'plane- wave' of infinite extent. 

There is no solution to the second problem if the 5-function in the amplitude 



(|3.6I) is interpreted as a Schwartz distribution. Consequently, the conventional 
solution to the third problem, i.e., the integration of the 5^-function in the dif- 
ferential cross-section (14.61) . which is explained in Sees. [5] and [61 and which is 
essentially equivalent to 'Fermi's golden rule,' is an ad hoc calculation unrelated 
to distribution theory whose physical justification is that it gives an answer that 



agrees with the data 

A mathematically consistent solution to all three problems is to redefine the 
transition amplitude (13.11) . and thus (13.61) . as a nonlinear generalized function in a 
Colombeau algebra. Instead of (|3.6I) we therefore write 



where (E) = p{E/e) /e is the scaled moUifier associated to the damper p in (17.11) 
which replaces the sharp cut-off introduced in Remark [H and where 1 /e corre- 
sponds to T/2 as remarked in Sec.|7l2. Therefore, the usual quantum mechanical 



rule (|4.1I) leading from amplitudes to probabilities yields 

where the limit T ^ oo, is taken at the end as in the conventional calculation. 

It remains then, just like in the conventional theory, to insure that the regular- 
ization implied by the damper p yields the proper physical result, i.e., Rutherford's 
formula. This has been done in Sec.[8]where, expressed as conditions on its inverse 



^As remarked in JT] p. 101], the use of 'wave-packets' to represent the incident and emerging 
particles implies that the appearance of squares of (5-functions is avoided. (See also, fTT, p. 105- 
106].) Mathematically, this is equivalent to chosing a smooth rather than sharp regularization of 
the Dirac ^-function, and working on the 'representatives' rather than on the 'distribution' itself. 
This bypasses problems 2 and 3, but does not solve them. 
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Fourier transform p, the conditions on p are given by (18.81) . i.e., 



p{z)eSiC), j dzp{z) = l, 
P(0) = ^' / 



J dz z"^p{z) = 0, Vn = l,...,g, 
dz |pp(z) = p(0). 



(9.3) 



The conditions on the first line of this set are necessary in order that the square of 
the 5-function whose representative is p2/T is well defined, so that it is possible 
to calculate the probability by squaring the amplitude. These conditions are 
equivalent to the statement that p2/T is an element of a Colombeau algebra Q, i.e., 
a nonlinear generalized function. 

The conditions on the second line are restricting the ^-function p, or equiv- 
alently its Fourier transforms p, to a particular subset of Q such that (19.21) yields 
the Rutherford formula. In combination with the conditions on the first line, they 
insure that the regularization implied by p corresponds to a special class of smooth 
wave-packets whose properties are equivalent to a plane-wave of infinite spatial 
and time extent 

As was remarked is SeclTJl, it is not easy to find a simple explicit example 
of a suitable moUifier p satisfying all conditions (|9.3I) . but there is not doubt 
that such moUifiers, and their corresponding dampers p, exist. The physically 
interesting question which is not addressed in the present paper is whether this 
class of moUifier/dampers is just a subset of Q such that first order time-dependent 
perturbation calculations make sense mathematically, or whether it corresponds 
to a physically significant redefinition of the notion of quantum-mechanical state 
with potentially far reaching implications for the foundations of quantum theory. 



^In the present paper the focus is on the time-integral in Eq. ( I3.3l l because the spatial Fourier 
transform of the Coulomb potential, Eq. ( I3.4l i. is well defined. In a general scattering problem, the 
integrations over all four dimensions are likely to need some kind of regularization. 
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